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Abstract 

In this work we study cat maps with many degrees of freedom. Classical cat maps are classified using 
the Cayley parametrization of symplectic matrices and the closely associated center and chord generating 
functions. Particular attention is dedicated to loxodromic behavior, which is a new feature of two-dimensional 
maps. The maps are then quantized using a recently developed Weyl representation on the torus and the 
general condition on the Floquet angles is derived for a particular map to be quantizable. The semiclassical 
approximation is exact, regardless of the dimensionality or of the nature of the fixed points. 
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1 Introduction 



Linearization of a dynamical system near a periodic orbit is one of the most fruitful starting points for the 
analysis of classical motion. In its turn, the symplectic group of linear Hamiltonian systems in plane phase 
space is easily quantized to form the corresponding metaplectic quantum group. Essentially, the generating 
function for the group of canonical transformations is simply exponentiated to obtain a representation of the 
quantum unitary transformation. 

If the chosen orbit is a point of equilibrium, the corresponding linear system belongs to the homogeneous 
symplectic group, characterized by a single equilibrium, usually taken as the origin. Likewise, the Poincare 
map in the neighborhood of a periodic orbit is linearized into a homogeneous symplectic map with discrete 
time. The essential character of the motion is classified according to the eigenvalues of the symplectic matrix 
M, that determines the evolution of phase space points x : 

x'^Mx. (1.1) 

There may be 

(a) pairs of eigenvalues (A, A^^); 

(b) pairs of eigenvalues on the unit circle (e*^,e~'^) 

(c) quartets of general complex eigenvalues A^^e**^. 

On varying parameters, it is possible to obtain unit eigenvalues , or eigenvalue collisions, but the above 
classification is generic for a given symplectic system |^] . 

It is always possible to decompose such a generic linearized system into sub-systems in invariant subspaces 
of two dimensions, corresponding to cases (a) and (b) above, or four dimensions in case (c). Case (b) is 
the elliptic map, which is trivially integrable, whereas case (a) defines hyperbolic motion. This is also very 
simple in the linear limit, but can become a source of chaotic mixing as nonlinear perturbations are added. 
Alternatively, this effect is achieved by wrapping the plane space itself into a torus. 

The resulting symplectomorphism of the torus is known colloquially as a cat map, characterized by a 
symplectic matrix with integer elements. A hyperbolic cat map is structurally stable, i.e. the orbit structure 
is invariant with respect to small nonlinear perturbations as a consequence of Anosov's theorem ||l|. The 
same is true of a four dimensional loxodromic cat map with general complex eigenvalues in case (c). These 
structurally invariant systems are known as Anosov systems, they are ergodic and mixing. 

It follows that four dimensions is the lower bound in which we can study loxodromic periodic orbits |^], 
characterized by stable and unstable manifolds where the orbits spiral inwards and outwards respectively, 
and their effect on the quantum energy spectrum. This is the reason for their absence in all previous 
studies of the quantization of cat maps, though Greenman has recently analyzed the periodic structure 
of higher dimensional classical cat maps. Dimension four is also the least dimension for the analysis of the 
decomposition of the neighborhood of orbits into elliptic and (real) hyperbolic components. 

In some cases this decomposition is only local, because the canonical transformation that achieves it is 
not itself a cat map. Then the quantum quasienergy spectrum will not be decomposed into the correspond- 
ing lower dimensional spectra. In any case, all cat maps derived from each other as a result of similarity 
transformations involving other cat maps are equivalent: they have the same (classical and quantum) eigen- 
values and the same number of fixed points. ( Note that, on the torus, a homogeneous linear map may have 
multiple fixed points.) 

For this reason, we discuss in the next section the integer subgroup of symplectic transformations, nick- 
named the feline group. For higher dimensions than two, we encounter the problem of a priori identification 
of a cat map. An alternative approach involves generating functions, of which there are several choices. 
However there is a great advantage to using generating functions that are invariant with respect to feline 
transformations . 

In section 3 we analyze the dynamics of classical cat maps and classify four-dimensional cat maps, pro- 
viding examples of various types. These examples are then quantized in section 4. They share a simplifying 
property that permits us to discuss the periodicity of the propagator and the exactness of the Gutzwiller 
trace formula, without analyzing the subtleties of general cat map quantization. This is the subject of section 
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5, where we determine the set of Floquet angles that allow the quantization of a given map; and study the 
feline invariance of the quantization. 



2 Feline-invariant generating functions 

A point in the even-dimensional phase space with L degrees of freedom on a 2i-torus has coordinates 
separated into L momenta and L positions, so that = ( ^ ) = ( '^'^ ). All the 2L coordinates 
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are periodic with periods Ag^ and Ap^. For simplicity we will treat the case where we can choose units so 
that A(7i and Ap^ are all equal to 1. The range of values of x is then the unit 2L-hypercube denoted from 
now on by □. 

Let us consider then a linear automorphism on the 2L-torus generated by the 2L x 2L matrix A1, that 

takes a point x- = i ^ ) to a point a:;+ = ( 

\ Q- J \ 1^ 

x+^Mx^ mod(l). (2.1) 

In other words, there exists an integer 2L-dimensional vector m = ( | , such that 

\mg J 

x+ = Mx^ - m. (2.2) 

The components of m denote the winding numbers made by the point X- around the respective irreducible 
circuit on the 2L-torus after the application of the map A4 . The torus will be divided into regions labeled 
by their respective vector m. For the map to be conservative, the Ai matrix must be symplectic, that is 



where is the transpose of Ai and 



M'ZM=Z, (2.3) 

(2.4) 






-1 " 


1— ' 






The matrix A4 must have integer coefficients for the 2L-torus to be mapped onto itself. 

For one degree of freedom ( L = 1 ) these systems are known as Arnold's cat maps [Q. If |tr(A^)| > 2 
the map has two distinct real eigenvectors, so it is ergodic, mixing and purely hyperbolic. For the case 
|tr(A^)| < 2 there are no real eigenvectors, the map is then elliptic. For |ir(A^)| = 2 there is only one real 
eigenvector and so the map is parabolic. For more degrees of freedom richer structure may appear, as will 
be discussed. 

The set of matrices A4 satisfying ( |2.3| ) form the symplectic group, so that the matrix A4' , obtained from 
M by the similarity transformation 

J\f-^MAf = M' (2.5) 

is also symplectic if the matrix M is itself symplectic and shares the same eigenvalues as M. Indeed, we may 
consider the symplectic transformation corresponding to A^' as the same as Ai, but viewed in an alternative 
symplectic coordinate system. 

Consider now a product of cat maps; this must be symplectic and all the matrix elements will be integers. 
Since the inverse of a cat map is also a symplectic integer transformation and the unit matrix likewise, the 
set of cat maps is a subgroup of the symplectic group, appropriately nicknamed the feline group. Again, we 



may consider similarity transformations of the form (2.5) among cat maps A4 and J\f as defining essentially 



maps viewed by alternative symplectic coordinates on the torus. 

The importance of defining equivalence classes of similar cat maps increases with the dimension of the 



torus. It is easy to define integer matrices and the a posteriori condition (2.3) merely restricts the value 
of the determinant to unity when L = 1. However, for L = 2, the symplectic property already implies ten 
independent conditions to be satisfied by the sixteen integer matrix elements. An alternative procedure, 
that we adopt here, is to define the transformation implicitly by means of a generating function, thus 
guaranteeing (|2.3|). The disadvantage is now that we need to ensure that the resulting matrix Ai has integer 
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entries, leading to additional restrictions on the allowed generating functions (see below ) . Nevertheless, 
this approach has proved very fruitful, because the quantization relies heavily on the generating function, 
and not on the transformation matrix itself. 

The generating function in position representation for the one degree of freedom case is 

S{q-,q+,m) = [M22ql - 2q^{q+ + nig) + Mii{q+ + rUgf - 2A^2iTOp9+] • (2.6) 

ZJVI21 

This function generates the dynamics through 

. dS{q-^1^^--) (2.7) 
oq+ 

= (2.8) 

oq- 

Thus the quadratic part of the generating function is common to the entire torus, whereas the linear part 
depends on the winding vector m that changes discontinuously on the boundary of each subregion of the 
torus. 

The weakness of the position generating function is that it transforms in a complicated way under feline 
similarity transformations. Only in the special case of a point transformation, that does not mix momenta 
with positions, will the generating function remain invariant. In contrast, the center and chord generating 
functions, are known to be symplectically invariant [Q. Adapted to the torus, we shall show that they can 
also be chosen to be invariant under feline transformations. 

The starting point is to define the center point x and the chord ^ such that 

x±=x±^^, (2.9) 



so that the center 

and the chord 

Given the initial point x- 



x+ + X- 



(2.10) 



2 

C = x+-x^. (2.11) 



X = liM + l)x_-^ (2.12) 

^ = {M-l)x.-in. (2.13) 

Elimination of x^, then establishes the direct relation between centers and chords: 

1{M + 1) 



(A^ + 1) 



e+(X-l)"'m (2.14) 
x-2{M + ir^ni. (2.15) 



Center and chord generating functions respectively denoted by S{x,ni) and S'(^,m) are defined in Q so 
that the transformation is obtained as 

» = (2.16) 

In analogy to the dynamics in the plane, we may consider that use of the center representation identifies 
the orbit with the reflection (or inversion) through x, because of ( 2.10| ). Hence we shall also refer to x as 
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the reflection point. The chord representation is locaUy equivalent to the uniform translation of phase space 
by the chord (2.11). Equating ( ^.14 ) with ( ^.16 ) and ( ^.15 ) with (2.17), we obtain the quadratic generating 
functions 



S{x,m) = xBx + x{B - :S)tn + f{m) 
where B and (3 are symmetric matrices; they are the Cayley parametrization of A4 : 



Jc/3C + ^^(/3 + a)m + 5(m) 



[l + M) 
[M + l) 
{M-l) 



(2.18) 
(2.19) 

(2.20) 
(2.21) 



If Ai has an eigenvalue equal to 1 then the /? matrix will be singular. This corresponds to a caustic 
of the center generating function. Whereas if M has an eigenvalue equal to —1, then B will be a singular 
matri x wh ich cor resp onds to a caustic of the chord generating function Q].Some useful relations obtained 
from ( |2.20| ) and (|^) are 



{B-Z) = -2Z{l+My 
(/3 + a) = -2^{M-iy 



(2.22) 
(2.23) 



1 



and 



M = 



(1 



ZB) _ {Z[3 + 1) 



(2.24) 



(2.25) 



(i+as) (a/3-1)' 

The functions /(m) and g(ni) are arbitrary, since they only depend on the winding number m, so they 
do not affect the transformation (2.14) or (^.15| ). However, the center generating function for cat maps can 
also be obtained directly from the map ( ^.2| ). This is a composition of the symplectic map on the plane , 
whose center generating function is Si{x) — xBx, with the uniform translation T-m of vector — m that pulls 
back the final point to the unit cell □. The generating function of such a translation is Q 5*2(2;) — ~mAx , 
also symplectically invariant. Then, using the composition law for center generating functions [Q we obtain 
the generating function (2.18) with 

/(m) = imSm. (2.26) 

As in the plane case generating functions are related among themselves by Legendre transformations. 
Thus, S{x, m) is obtained from the more familiar position generating function S{q-, (7+, m) as 



S{x,in) = 5((?_,g+,m) + -(p_ +p+){q+ -q^), 

whereas the relation between the chord and center generating function is 

S'($,m) ^iAx-S{x,m). 



(2.27) 



(2.28) 



In each case, the variable absent on the left s ide is eliminated by requiring the right side to be stationary 
with respect to it. The skew product in (2.28) , 



(2.29) 



is the symplectic area of the parallelogram formed by any pair of vectors ^ and rj. Then, using (2.28), for 
the center function with the term ( 2.26| ), we obtain the chord generating function ( ^.19 ) with 



g{m) — -m/9m. 



(2.30) 
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We shall label periodic points of period I as xi. Thus, fixed points xi are such that the chord ^ = 0, or 
the center x ^ xi = {A4 — l)^^m, which inserted in (2.28) leads to 



= 0,m) -S{xi,m). 

There follows the restriction that the terms in 5(^, m) and S{x, m) that depend only on m satisfy 



(2.31) 



/(m) + .g(m) = -ni(/3 + B)m.. 



The choice ( |2.26 ) and (2.30) obviously satisfy this criterion, but another possibility is 



/(m) = im(B + 3)m and ^(m) = |m(/3 - 3)m, 



where we define the symmetric matrix 



" 


1 " 


1 






Using (2.33), we match the value of the action for a fixed point previously proposed by Keating 
position representation 

S{q- = qf,q+ = qf,m) = S'(a;i,m). 
In conclusion, the center and chord generating functions for multidimensional cat maps are 

S{x,m.) — xBa; + a;(_B — 3)m + ^m(_B + -3)m 
5(^,m) = l^/3e + i^(/3 + a)m+im(/3-5)m. 



(2.32) 

(2.33) 

(2.34) 
for the 
(2.35) 

(2.36) 
(2.37) 



The corresponding generating functions for the transformation = A4x- in the plane are just S{x,0) 
and S{^,0). It is important to note that all the reflections of the torus can be obtained with center points 
whose coordinates are in [0, ^]. It would thus be possible to define such center points x' as 



+ X- 



mod(-). 



(2.38) 



This choice would not lead to the explicit relation (2.12) with the winding number of the transformation, 
but instead 

x'^^iM + l)x^-^, (2.39) 

where m' has coordinates — (m^ or m.i+1), so that all the coordinates of x' would be in [0, Hence, we 
allow instead the center points x, defined as in (2.10), to have coordinates in the full interval [0, 1], keeping 
the explicit relation with the winding number of the transformation. 

As a consequence, the chords defined as in (2.11), have coordinates lying in the extended range [—1, 1]. 
Thus, center points differing by integer loops around the torus are equivalent and so are chords differing by 
two integer loops. 



X = x + k 

e = e + 2k. 



(2.40) 
(2.41) 



We find that ( 2.40| ) and ( 2.41 ) imply that, in ( 2.14 ) and ( 2.15 ) respectively, the winding number m is 
equivalent to: 

m = m' = m + (7W - l)k (2.42) 



in ( |2.14| ) and 



m = m" = m - (7W + l)k. 



(2.43) 
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in ( ^.15 ). This implies that replacing m by m' in the generating function 5(^,m) will lead to equivalent 
center points related by (2.40). To obtain equivalent chords related by (2.41), it is necessary to replace m 
by m" in the center generating function S{x, m). Performing the mentioned replacements we will obtain: 



5(^,m') 
S{x, m") 



S{^, m) + e A k-imrik-ikAik 



1 



1. 



S{x,-m) - 2a; A k--mr2k--kA2k, 



where 



Ai 

r2 

A2 



+ {Z + Z]M 



(2.44) 
(2.45) 

(2.46) 
(2.47) 
(2.48) 
(2.49) 



In this way we can restrict m to integer component vectors that lie in one of the two fundamental 
parallelepipeds 



{M 



!)□ 
-!)□ 



for^(e,m) 
for S{x, m) 



(2.50) 
(2.51) 



where □ is the unit hypercube that denotes the 2L-torus. Hence, the different orbits denoted by a given 
chord ^ are given by all the integer m lying in Oj. The number of such orbits is independent of ^, so taking 
^ = 0, we equate this to the number of fixed points tj, i.e. 



TiM) = |det(X - 1) 



-,2L 



|det(3/3-l)| 



(2.52) 



The different orbits that have the point x as its center are denoted by all the integers m lying now in Oj 
The number of these orbits is given by the volume of O^; which is 



t{-M) = |dct(M + l)| = 



■)2L 



det(aB + l) 



(2.53) 



Note that the number of periodic points of period two is t^t^;. 

For the matrix A4 to represent a 2L-cat map it must be symplectic, so that the map is area preserving, 
and A4 must have integer entries. Let us now translate both the conditions for the B and (3 matrices. The 
first condition implies that B and f3 are symmetric matrices, and any symmetric matrix is associated through 
( 2.25 ) to a symple ctic m atri x. Th e second condition restricts the B and /3 matrices to have rational entries. 
Indeed, following (|2.20D and(|2.2lD we will have 



B = 



P = 



B _ B 
det{M + 1) ^ 

det{M - 1) 



(2.54) 
(2.55) 



where B and (3 are symmetric matrices with integer entries and the denominators are defined by (2.52) and 
(2.53). It can happen that all the coefficients of the matrix B { or (3 ) have a common factor that is not 
coprime with (respectively ). Dividing by this common factor the fraction in (2.54) ( or in ( p. 55 ) ) is 
reduced to the form 



B 



B 



(2.56) 
(2.57) 
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But not any symmetric B or (3 matrix with rational entries guarantees that the associated symplectic matrix 
will have integer elements. So we must find conditions on B and /? for this to occur. 

The characterization of the matrices B ov (3 requires that the corresponding transformation on the plane 
maps the points of an integer lattice among themselves. We will examine the case L = 1, as the extension 
for many number of degrees of freedom follows easily. There are two fundamental chords corresponding to 
fixed points on the torus: 

6 = ( J ) and 6 = ( J ) , (2.58) 

leading to the fixed points: 

Xj = 1(3/3 + ^ with j = 1, 2. (2.59) 

Of course, there is also = 0, but this "plane fixed point" makes no restriction on the torus map. For the 
transformation to be a cat map, all the corners of the fundamental parallelogram □ must be fixed points, 
so, for any integers r and s, there are integers mi and rn2 such that: 

mia:i + 1712X2 = ( I ) = ^(3/5 + 1) ( ) • (2.60) 

This is true if only if 2(3/3 + 1)^^ has integer entries. This condition is general for any degrees of freedom. 
Similarly, we find that 2(3/3 — 1)~^ having integer entries is also a necessary and sufficient condition for the 
corresponding to define a cat map. But , if M defines a cat map, so does —M with the associated 

chord matrix —B. Therefore, it is also a necessary and sufficient condition, for a center generating function 
to determine a cat map, that the associated center matrix B have the property that 2{ZB ± 1)~^ be an 
integer matrix. Evidently we easily find a subclass of cat maps by restricting their Cayley parametrization 



to B (or /3 ) matrices of the form (2.54) such that det{l ± ^B) = ±1,±2, then 2(1 + ^B)^^ is an integer 
matrix. 

Although the conditions on the symmetric matrices i? or /3 to denote a cat map are not as trivial as the 
ones on the symplectic matrix , it is simpler to find rational symmetric matrices that fulfill the condition 
on /3 and B than to find integer symplectic matrices. The fact that B ov (3 are symmetric and of the form 



( 2.54 ) allows us to find cat maps by sampling [(L) x (2L + 1) + 1] integer numbers. Otherwise, to fulfill the 
condition (2.2), needs a loop over (2L)^ integer coefficients. 



To conclude this section, we verify the property of feline invariance for the chord and center generating 
functions. First, we note that, symplectic invariance in the plane jj] implies that under a symplectic coor- 
dinate transformation x ^ x' ^ J\fx, S{x,0) = S{x',0) and S{^,0) = S'(f',0), with = Af^. But it is also 
evident that the winding number m transforms in the same manner: m' = Mm. As far as the x dependent 
term in S(x,m), we thus find that the effect of the feline transformation is merely that of substituting 
B — > J\f*BJ\f, and similarly the change in 5'(^,m) is obtained from /3 M*(3N. The constant terms /(m) 



and g(m) in (2.18) and (2.19|) are not invariant under a feline transformation in the form (2.33) that we have 



chosen to match reference p], so that it is preferable to use (2.26) and (2.3C) when dealing with equivalence 
classes of cat maps. 

It is important to note that, unlike the symmetric matrices B and /3, ZB N~^ZBN and 3/3 
Af~^^PAf, under a similarity transformation Ai Af~^AiM. Therefor, the eigenvalues of ZB and 3/3 are 
feline invariant, just as those of A^, and can thus be used to classify cat maps. 



3 Classification of classical cat maps 

The periodic orbits of cat maps have been studied in great details by Pcrcival and Vivaldi |^ and also by 
Keating in for one degree of freedom and their results were recently extended to an arbitrary number 
of degrees of freedom It is shown that a point on the unit 2L-torus is periodic if and only if all its 
coordinates are rational and any grid of points with rational coordinates is invariant under the action of 



the map. From (2.2) we can see that the periodic points xi of integer period I are labeled by the winding 
numbers m, so that 

'M=(M'-l)-im^3/3W-l)f. (3.61) 



xi 
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Here fi^''^ denotes the symmetric matrix associated to through ( ^.21 ). To have xi on the unit 2L- 
hypercube □, m must he within the parallelepiped formed by the action of the matrix (A^' — 1) on □. Hence, 
the number of integer points m is given by its hypervolume, so that the number of periodic points with 
period / is 

22L 

.(A.') = |det(A.'-l)H^j^^^^^^^. (3.62) 

According to ( 3.6l| ) the periodic points of period I form a lattice in phase space with rational coordinates. 
The motion of any point x- = xi +5- near a fixed point xi will be 

Mx^ = M{xi + =xi+ M5- =xi+5+= x+. (3.63) 

To determine the character of such a motion we have to study the eigenvalues and eigenvectors of the matrix 
M: 

>^4 = \>^iW'' (3.64) 

The modulus | | indicates that the motion is stretching while the argument e'^'' indicates rotation around 
the fixed point xi. For a symplectic matrix M, if Xjvi is an eigenvalue of M then A^, y— and -p— will also 

be eigenvalues of M. 

The classification of the eigenvalues is possible in either the A4,B, or f3 descriptions. Using (2.25) we 
obtain the relation with the eigenvalues of and -3/3 denoted respectively as A^js and A3/3, 

(l-Ajg) (A3^ + 1) 

= (TTa^ - (A^ ^'-''^ 

and inversely: 

(1-Aa^) fo 

Aas = — — — - (3.66) 

(i + Am) 

> i^M + 1) 1 

[Am — tj AjB 

In this way, if A^b is an eigenvalue of ^B, then A^^, — Ajs and — A^^ will also be. In the same way, if A^/3 
is an eigenvalue of 3(3, then AJ^, — Aj/3 and — A^^ also are. 

For cat maps with L = 2 the matrices M, B, and /3 will be 4 x 4 and we then have the following generic 
cases for the eigenvalues 

1. Elliptic: there is a conjugate pair of Xm both on the unit circle and conjugate pairs of purely imaginary 
X^B and A3/3- 

2. Hyperbolic: there is a pair (A^^, j^) on the real axis and pairs {X^b, ^A^b) and (A3/3, — A^^) on the 
real axis. 

3. Parabolic: there are degenerate eigenvalues Xm — ±1- Then for Xm = 1, /3 is singular and Ajs — 0; 
for Xm = —1, B is singular and A^^ — 0. 

4. Mixed: each pair belongs to a different one of the above categories. 

5. Loxodromic: the eigenvalues of A^, _B. or /3 matrices form quartets of complex eigenvectors of the form 
(^■'^^ AXi"' '^A4' A^) ^^"^ matrix and (A^s, —X^s, AJg, — A^g) for S being one of the symmetric 
matrices B or /3. 

The first three cases arise also for one degree of freedom cat maps. The case L = 2 is the lowest number 
of degrees of freedom where not only elliptic, parabolic and hyperbolic fixed points will appear, but also 
loxodromic ones. For more degrees of freedom no new cases will occur; there is only a greater variety of 
mixed cases. Nongeneric possibilities arise for any dimension for continuous families of systems ||l[, but we 
do not know if there exists any corresponding cat map. 
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In the general case, the motion is the composition of stretching and rotation, but if the angles of the 
rotation are of the form 6k = j27r, after j applications of the map, the matrix M.^ has only real eigenvalues 

that come in pairs |A^|"' and |A^| "'. The dynamics has then an ignorable coordinate; the angles are 
constants of the motion for Ai^ . For one degree of freedom systems, where there are only two eigenvalues, 
\m and A^, the condition 

Tr {M) = \m+>^*m= integer (3.68) 

implies that in the elliptic case only angles of the form Ok = jtt with j = 2 or 3 are allowed. This is an 
example of a rational restriction on 9k, that is, Ai^ will be the identity map. This result is in accordance with 
Marie's theorem that two-dimensional symplcctomorphisms are either Anosov or they have zero entropy. 
Nonetheless, irrational rotation angles do exist in the loxodromic or mixed cases for i > 1. 

We now turn our attention to the classification of four-dimcnsional cat maps. Recalling that the charac- 
teristic polynomial for any k x k matrix A is 



P^(A) = det(A-A) = ^a„A* 

n=l 

where the q;„ coefficients are given by the recurrence relation, 

ao = 1, 



1 



an-iCii with tti = Tr(A^), 



n ^-^ 

i=l 

we obtain that for any symmetric matrix B (i.e. B and /3 ), 
so that 



Pas(A) - A^- -&2A2 



detS 



(3.69) 

(3.70) 
(3.71) 

(3.72) 
(3.73) 



with 62 = Tr 



(dB) , ( note that dct^B = det B). There is a similar expression for P3/3(A), required when 



the center representation is singular, 



Po;3(A) = A^- -/?2A2+det/3. 



(3.74) 



where (32 = Tr 



For the symplectic case, we obtain 



Pm{>^) = A'* - Tr{M)\ + - [Tr{M^) - Tr{Mf] A^ - Tr{M)\^ + 1 



(3.75) 



which is harder to analyze, so we will perform the classification of the different behaviors using or 3/3- 
Solving P3_b(A) = using (3.73), leads to 



\ 



whereas 



he = ± 



\ 



(3.76) 



(3.77) 



We can now classify the different behaviors according to the feline invariants of the matrix ZB or 5/3, the 
eigenvalues of the symplectic matrix M. being obtained with the help of ( 3.65| ). The loxodromic behavior 
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corresponding to four complex eigenvalues will appear if the square root term 
square root is negative. That is, the loxodromic behavior will appear only if 



If) -deti? 



detB > ( ^ 



or equivalently 



det/3> 



inside the 



(3.78) 



In fig 3.1 we find a complete classification of the different types of behavior according to the feline 
invariants of -3/3 and the same arises for the invariants of ^B. These invariants also allow us to obtain r^, the 
number of orbits for any chord ^, 



or Tx the number of orbits centered on x as 



24 




2-4 det B 












24 det p 


^35(1) 







(3.79) 



(3.80) 



Examples 

We here show some examples of different types of cat maps. The examples are such that the (3 matrix, that 
denotes the chord representation, has integer elements. 

The Hannay and Berry Cat 

It is important to see that the first cat map to be quantized ( the Hannay and Berry cat map) [Q can also 
be treated with the formalism described here . For a 2 x 2 matrix, the characteristic polynomial reads, 



For the symplectic matrix 
Mhb -- 

so that 



2 1 

3 2 



with the associated symmetric matrix (3hb 



3 
-1 



(3.81) 



(3.82) 



(32= Tr 



the number of fixed points is 



hb 



iZPhb) 

22 



PwW 



-6, 



= 2 



hb 



= 2 + V3 
= 2- V3. 



and the eigenvalues of the matrix Mhb are 

The map is then hyperbolic. 
The double hyperbolic case 

Let us now study cat maps that have two pairs of eigenvalues both in the real axis, that is, maps belonging to 
the first quadrant in figure 3.1. The (3hh matrix below describes this case, the associated symplectic matrix 
being Mhh ■ 



Phh — 



2 12 
2 2 1 
12 
2 10 1 



and M hh 



2 -2-10 
-2310 
-12 2 1 

2-201 



(3.83) 
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MIXED 

Figure 3.1: A classification of the different cat maps for two degrees of freedom according to the feline 

invariants of the ZP matrix. The parabolic boundary is det/3 = • A similar picture exist for the ZB 

matrix. 
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The invariants are 



det Phh = 17 and /32 = Tr {^fihh) 



= 20. 



Thus, the number of fixed points ~ 2 and the eigenvalues of the symplectic matrix M are 



Af' = 2.112388 



A 



hh 
2 
\ hh 
^3 



A 



hh 



0.4739 

5.22 

0.1914. 



The mixed case 

We consider here the case where the eigenvalues of the map are of mixed nature, i.e. there is a couple of real 
eigenvalues and a complex conjugated pair on the unit circle. This corresponds to /3 matrices that belong to 



the third and fourth quadrant in figure 3.1. 

One matrix of this kind is Pehi, the associated symplectic matrix being Aiehi 



Pehl — 



and Mehi = 








-1 

2 
1 





-1 
1 
2 



(3.84) 



The invariants are 



detPehi - 

Thus, the number of fixed points t^^^ 



-15 and 132 = Tr 



{ZPehlY 



4. 



1 and the eigenvalues of the symplectic matrix M. are 



Afi 

\ehl 



\ehl 



A 



'.hi 



= 2.6180 
= 0.381966 

= expK- 



27r 

exp 1 -t — 



This example shows rotation angles that are fractional multiples of tt, for which the dynamics are equivalent 
to that of hyperbolic systems with one degree of freedom . Another example of a mixed system is given by 
the matrix Peh2, the associated symplectic matrix being Aieh2 ■ 



eh2 



and M, 



:/l2 





-1 


2 



~1 
-1 

2 

2 





-2 
1 
3 



Now we have 



= 4. 



det/3e/i2 = -7 and /32 = Tr {ZPeh2) 
Thus, the number of fixed points r|''^ = 2 and the eigenvalues of the symplectic matrix A4 are 



(3.85) 



(3.86) 



A? 



eh2 



\eh2 

^2 

\eh2 



A 



;/i2 



3.0906578 

0.32355571 

exp (a. 27354496) 

exp (-il. 27354496). 



In this example there are no rotations with angles that are a fraction of tt, the dynamics will then be ergodic 
and mixing in the whole phase space. 
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The loxodromic case 



We choose now a (3 matrix that belongs to the loxodromic region in Fig 3A, for example (3io 
associated symplectic matrix being Mehi ■ 



lox 



2 10 
2 2 1 
12 11 
111 



and M. 



lox 



10 

110 

1-111 

-1 -1 -2 



whose 



(3.87) 



whose invariants are 



det/3/o2: = 5 and (32 ^Tr {ZPiox) 



= -4. 



(3.88) 



Thus, the number of fixed points r^°^ = 2 and the eigenvalues of the symplectic matrix M. are 



\lox 



^3 
\ lox 



1.7000157exp(il. 1185178) 
1.7000157exp(-il. 1185178) 
0.5882298 exp(il. 1185178) 
0.5882298 exp (-il. 1185178) . 



There are no rational rotations, so the dynamics is ergodic and mixing in the whole phase space. 



The double elHptic case 

In this case, the eigenvalues are characterized by two rotation angles. An interesting feature is that the 
condition that the M. matrix have integer elements, prevents either angle from being an irrational multiple 
of TT. This restricts the dynamics to the trivial integrability observed for one degree of freedom systems. 
In conclusion, any symplectomorphism on the 4-torus belongs to one of the following cases: 

1. Ergodic on the full phase space (a manifold of 2 degrees of freedom). For the double hyperbolic, mixed 
or loxodromic cases with irrational rotation angles. 

2. Ergodic in a 2-dimensional manifold . We have a constant of the motion, then the ergodicity is 
restricted to a low dimensional manifold (a manifold with L ~ \). For the mixed or loxodromic cases 
with rational rotation angles. 

3. A root of unity: that is, we have two constants of the motion. All orbits are periodic and have the 
same period for this double elliptic cases. 



4 Simple Quantum Cat maps 

Quantum dynamics is characterized by a unitary evolution operator, or propagator. It is possible to obtain 
the center and chord representation of an operator on the torus from its counterpart on the plane. This will 
be the way that we quantize cat maps, leaving the general construction for the next section. In some cases 
the quantum propagator denoting cat maps in the center or chord representation acquires simple expressions; 
in this section we will treat these special cases. More details about torus Hilbert space and its center and 
chord representations based on reference are available in Appendix A. 

It is well known that torus quantization implies that the Hilbert space [tC^]^ associated to the 2L-torus 
has finite dimension and is characterized by a vector Floquet parameter x — {Xp^ Xq) whose components 
are real numbers belonging to [0, 1] . The fact that has finite dimension, implies that position and 

momentum eigenstates can only take on a set of discrete values that form a discrete lattice called quantum 
phase space (QPS). Any point x in this QPS has coordinates, 

x=( ^™ \ =L(^ + Xp\^ (4.1) 
\ q-a ) N \ n + Xq ) ^ ' 
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The center and chord representations are based on translation and reflection operators on this QPS as 
we show in [p| • Chords arc of the form 



with r and s integer numbers and ^ — 
a and b, 



1 

N 



1 - 



(4.2) 



while the center points Xa,b are labeled by half-integer numbers 



Xa,b - ^ 



a + Xp 
b + Xq 



(4.3) 



At a first stage we restrict our attention to those maps that can be quantized on the Floquet parameters 
X — (0, 0). As we will see in section 5 this implies that the matrix Ai must satisfy 



''^MijMij+L = even integer for all 



(4.4) 



As we will show in section 5, if 2r^, defined in ( 2.57 ), and N are coprime numbers, a complete represen- 
tation of the propagator is obtained having the symbol on a lattice of chords S such that 



N ' 



(4.5) 



where the components of ^ are integer numbers up to N. So for any chord ^ there is an equivalent chord ^. 
For the allowed values of N, the propagator for cat maps in the chord representation takes the form 



U^(S) = 2^ (re) 



-i27rAr[iH/3E 



-i27rArim(/3-3)r 



(4.6) 



where is the number of fixed points of the classical map. For f3 m atrices that fulfill the feline conditions, 
the symbol U^(S) must represent an unitary operator. In that case ( A. 39 ) shows that it must has the form 



U^(S) 



~^27^N[^S|3= 

T ^ 



N 



which restricts 



N 



(4.7) 



(4.8) 



meOf 



At first sight, the phase cpNi-M) is only an unimportant global phase factor, but the interference of the 
different (^7v(A^0 for the different powers I of the map will have a crucial importance for the density of 
states. 

As ^ and S are equivalent chords, the symbol U»(^) and U»(S) are related by symmetry relations 
( |A.3l| ) so that 



Ua,(0 



-i27TN[S(C,n)] 



N 



(4.9) 



where S{^, n) is the action of the classical orbit whose chord is ^ and that performs n loops around the torus 
as defined in ( 2.37 ). 

If we chose the symplectically invariant form (2.30) for the ^-independent part of the chord generating 
function, instead of (2.33), we would have in (4.9) a supplementary phase factor g»2TArin0n _ gJ7„ ^[^]^ 
a "Maslov index" for the orbit. This observation is true for all the following quantum theory. 

In the case of the center representation, for 2t^ , defined in (2.56), and N coprime numbers, a complete 
representation of the propagator is obtained by performing a transformation to center points X that are 
integer multiples of 

X^x+^i^^X, (4.10) 
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where the components of X are integer numbers up to Q . On these points the center representation of 
the propagator takes the form 



i^t27TN[XBX\ ^ g 



r 



^iip'^(M)^2TTN[XBX\ 



(4.11) 
(4.12) 



where the last equahty is obtained by imposing the unitarity of U^vi and using (A. 40). Hence, we define the 
angle (p'j^{M) so 

^iip'^(M) ^ 2^Tx^ ^ ^2TTN\in(B+Z)m_ (4-13) 



From the symmetry relations ( A. 34 ), we find that the symbols on the original points x are 



(4.14) 



where here S{x,]) is the center generating function, defined in (2.36), on a center point x for an orbit 
performing j loops. The cases above are then special cases where the propagator on the torus has the same 
form as its equivalent on the plane; thus, they are ideally suited for the comparison of classical and quantum 
motion. 

To obt ain the more familiar position representation of the propagator from its chord representation, we 
use (|A.35|) , 



UAi(q„,qr! 



^ exp \ -i2nN 



n) + 



While ( A.36| ) allows us to obtain the position representation from the center one 



(4.15) 



U7w(q™,qn) = — — 2^ exp|i27riV p(Xp , j)+ (q,„ - g„) XpJ | 



(4.16) 



For the case of one degree of freedom, this leads to the Hannay and Berry propagator. 

As we will see in section 5, condition (4.4) is preserved for the different powers of the map, hence, if a map 
is quantizable on the Floquet parameters x = (0; 0)i all its powers also are. Furthermore, the propagator for 
A^' has the form 

UME)^'- ^e-2-^[^^^"'^], (4.17) 



N 



in the chord representation, where denotes the symmetric matrix associated to A^' through (2.21) and 
S are points on a lattice of the form 

„_r(MO^ 



N 



(4.18) 



As we discussed in section 2, lattices of rational points are invariant under classical cat maps. The 
denominator g of these rational points can then be used to label the lattice, i.e. there exists a minimal 
period Ig under which all the points on the lattice are fixed for A^'s. In other words, Ai''^ , restricted to the 
mentioned lattice, is the identity. We will call Ig the classical periodicity function of the lattice. As we can 
see in the appendix, torus quantization is performed on a lattice of rational points whose denominator is N. 
Under these considerations we find that the quantum propagator is also periodic. 

It is now possible to define conditions on the matrix Ai and the dimension of the Hilbert space N under 
which Um reduces to the identity operator, i.e. 



U 



M 



Are ^. 



This is best seen in the center representation, where we have (A.38) 



(4.19) 
(4.20) 
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with ,fN{x) defined in (A. 28). For the case where N is an odd integer, we can perform the quantization on 
points X that are integer multiples of ^ , as we already discussed. We can see that if the B matrix has 
all its elements that are multiples of p-, then the propagator given by ( 4.11 ) has the form ( 4.20 ), where 
fN(X) = 1 and (p — Lp'pf{ M). H ence, the B matrix denotes the identity if all its coefficients are multiples of 
This implies through ( ^.25 ) that the A4 matrix must have the form 



M 



mod(iV), 



(4.21) 



in agreement with reference [Q for the L = 1 case. 

For any matrix Ai and for a given odd value of N, there is some smallest integer k{N) such that 



= 1 mod{N) 



and in this way 



TJ" 



k{N) 



le 



(4.22) 



(4.23) 



In this sense, we may call the quantum propagator periodic with a period equal to k{N) that we then c all th e 
quantum period function (QPF) although it is completely determined by the classical map through (4.22). 
Note that k{N) = Zjv, i.e., the quantum and classical period function coincide for N an odd number, as we 
would expect. It is now easy to see that the eigenangles dm of the unitary propagator XJm are then 
restricted to lie on the k{N) possible sites 



2jTr + (p{N) 
k{N) 



,l<j< k{N). 



The eigenangle spectrum may be related to the traces of the propagator 



Tr 



m=l 
k{N) 



(4.24) 

(4.25) 
(4.26) 



where dj is the degeneracy of the jth site defined by (4.24). The density of states, 

i—l 1 — — 00 

is clearly invariant under 9—^9 + 2t:. Use the Poisson summation formula, leads to the trace formula, 

27r TT 



1 °° 
P{0) = ^ + -RcY^Tr 



1=1 



U 



M 



(4.27) 



(4.28) 



which holds for all maps on the 2L-torus. Now, for the cat map, the fact that the propagator has periodicity 
k{N) means that the density of states can also be written in the form [O 



k{N) 



(4.29) 



1=1 



That is, the eigenangles are restricted to the sites aj in ( 4.24 ), with the degeneracy at the jth site being 
given by 



k{N) 



(4.30) 
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This last equation then leads to simple expressions for several important properties of the eigenangle distri- 
bution; for example 

k{N) , k{N) r n 2 



(4.31) 



It is important to note that the (5-functions appear explicitly in (4.29), because the propagator is periodic. 
The different powers of the quantum map in ( |4.29| ) contribute only to the finite degeneracies, not to the 
(5-function form of the trace formula. Moreover, for all finite N, only a finite number of these powers are 
required to determine the degeneracy at every available site, and hence give the whole spectrum. However, 
k{N) —^ooa^sN—^oo, thus more terms are required to determine the spectrum as the semiclassical limit is 
approached. 



We now take the trace of ( 4.17 ), using ( A. 37 ), to obtain 



Tr 



UM.(e = o) 

2^(r(X'))-^ J2 exp 



-i27r7V-m(/3(') - 3)m 



where we recognize the action of the periodic orbits in the summation exponent, so 

n 



N 



(4.32) 
(4.33) 



Tr 



2^|det(^'-l)f * J2 e-^2^^KM 



N 



(4.34) 
(4.35) 



Note that m g implies that we are summing over the periodic orbits of period I. Equation (4.34) 
represents the Gutzwiller- Tabor trace formula [ p^ for the cat map, though we must note that, instead of 
being a semiclassical approximation, it is exact in this case. We r nust also note that the expression ( 4.35| ) 
for the traces of the propagator implies that in the expression ( 1.2£ ) for the density of states, we have 
interference of the different phase factors e*'''"*^-^'-' for the different powers I of the map. 



Examples : 

We will now study some quantum features of the classical examples studied in section 3. All of them fulfill 
condition (4.4), so that their symmetric /3 matrix has integer entries. Then quantization can be performed 
for all odd values of N, for which the quantum propagator will have the form 



U^(S) = ,/^e— [i-'^-] 



(4.36) 



j^. In the 



in the chord representation. In this case, the chords S, defined in ( 5.70| ), form a lattice of spacing 
following, we will study the quantum period function for these maps. 

We take cat maps of two degrees of freedom already studied in sectio n 2. These include all the possible 
classical behaviors and we will study their effect on the QPF defined in (4.22). 



The Hannay and Berry cat map 

To compare the different types of behavior that occur for two degrees of freedom, it is important to present 
t he a lready known QPF for the Ha nna y and Berry cat map whose symplectic matrix is Adht defined in 
( 3.82 ). The QPF is shown in figure 12, where we can then see that, although it has an irregular behavior 
as a whole, most of the points lie in families of straight lines which admit a maximum slope. Indeed, this 
indicates that there exists an increasing sequence of primes p such that < C for some C independent 
of N. According to Degli Esposti et al this implies that the quantum map is ergodic and mixing in 
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the semiclassical limit, within the definition of Von Neumann The average behavior of the QPF was 

estabhshed by Keating ll^ and it is shown that in the semiclassical limit {N oo) 



^ N / 



(log log N) (log log log N) 



(4.37) 



that is, the degeneracy grows very slowly with N. 




Figure 4.2: The QPF for the Hann ay a nd Berry cat map, chaotic with one degree of freedom, whose 
symplectic matrix is Aihb, defined in (3.82). 



The double hyperbolic case 

We now choose the symmetric matrix f3hh whose associated symplectic matrix is given by A4hh defined in 
(3.83). The QPF of this map is shown in figure 4.3. As we can see, the behavior is very different of that 
obtained for the Hannay and Berry Cat map. There are now families of parabolas instead of straight lines. 
The role played by N in figure 4.2 is here played by iV^, because there are here N'^ states for this system. 



We conjecture that this kind of behavior implies quantal ergodicity and mixing for two degrees of freedom 
systems in the semiclassical limit, although a more formal study of this conjecture will be realized in a future 
work. 



The mixed case 



We first consider the map whose symmetric matrix is Pehi, defined in (3.84), whose QPF is shown in figure 
4.4. This is compared to /3eh2, defined in (3.85), characterized by irrational rotation angles, whose QPF is 
now shown in figure 4.5. There is a marked difference in the behavior between figure 4.4. and figure [f.5| . 
While the former is very similar to the one obtained in figure 4.2 for a chaotic one degree of freedom system, 
the one of figure 4.5. is close to figure 4.3. for a system of two degrees of freedom. The very different 



behavior, as was explained in section 3, is due to the fact that the classical eigenvalues Xf^^ and Xf^^ denote 
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Figure 4.3: The QPF for the double hyperbohc map, whose symplectic matrix is Mhh, defined in (3.83) 



rotation by angles that are a fraction of tt. Then the behavior of the system will be equivalent a hyperbolic 
system with a single degree of freedom and it is ergodic only in a subspace. On the contrary, and Xf^^ 
denote rotation angles that are irrational fractions of tt, so that classical and semiclassically the ergodicity 
appears in the whole phase space. 
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Figure 4.4: The QPF for the mixed map , with rational rotations, whose symplectic matrix is A4e/ii, defined 



in (3.84) 



15 s^. 



2 O O O 3 O O O 



5 O O O 



Figure 4.5: The QPF for the mixed map, with irrational rotations, whose symplectic matrix is M.eh2^ 



defined in (3.85) 
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The loxodromic case 



Let us no w st udy the map denoted by the Piox matrix defined by ( 3.87 ). The corresponding QPF is shown 
in figure 4.6. The behavior is then similar to the one obtained for ergodic systems with two degrees of 
freedom. This map then manifests classical and semiclassical ergodicity on the whole phase space or the 
corresponding Hilbert space. Similarly to the previous example, we found that rotation angles that are a 



fraction of tt behave similarly to hyperbolic one degree of freedom maps shown in figure 4.5 



^ o o o 



Figure 4.6: The QPF for the double loxodromic map of two degrees of freedom, with irrational rotations, 
whose symplectic matrix is Miox, defined in (3.87). 
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5 The quantum feline group 



In this section we present the general quantization of multidimensional cat maps constructed in section 
two. In a first step we establish the conditions for the maps to be quantizable and then we construct their 
center and chord representations, based on the formalism developed in reference j^] and summarized in the 
Appendix. 

The classical automorphism generated by A4 in plane phase space being linear, its quantization on the 
Hilbert space [TYr]'^ , associated to the Euclidean phase space, will have the crucial property that 

UmT^uI, = (5.38) 

for any translation operator Tj. We now show how to associate to any Um ^ unitary operator on [7i^]^, the 
torus Hilbert space characterized by the Floquet parameter x- For this purpose let us study the restriction 
of Um to [H^]^ . We then have: 
Proposition 



Um [n^f C [H^\ (5.39) 
where 

N 

X' Mx- Y^iM^M) mod(l) (5.40) 
where we have defined the vector product of two matrices as the vector whose components are 

L 

{A^B)^^ = Y,A,B^,J+L (5.41) 



L 



Proof: Equation (5.38) implies 

UmT^ = Tm^Um- (5.42) 

We now restrict considerations to chords that perform integer loops around the torus, ^ = m, with m an 
integer vector with 2L components, i.e., chords that generate translations that are equivalent to the identity 
on the torus. We thus act on a given state of [H^]'^ on each side of (5.42), so that inserting (A. 20), we 
obtain 

for any vector m. Now we chose n = Aim , hence, we find that 

nAx' = A4-inAx + (A^"i)*5(A^"^) -5 n mod(l) (5.44) 

and we specify 2L independent integer vectors n = Cj {j = 1, 2L), such that each of them denotes one loop 
around one of the j-irreducible circuit on the torus. In this case CjZcj — 0, using the symplectic property of 
Ai (2^), and the fact that ~ 5, we obtain 



ejhx! = M ^e-j hx+^e.jZMZM^Zej'OYodil). (5.45) 



Remarking that 

^ {ZMZM'Z)^ = -^Y^^-^i^W. (5.46) 



we finally obtain ( ^.40 ) after straightforward manipulations. 

In the case L — 1, ( ^.40| ) reduces to the already known result , ||l| 



2 ^ V -^^21X22 
For two degrees of freedom [L = 2) we obtain 



^-Mx-^Z( ] mod(l). (5.47) 
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N 



mod(l). 



(5. 



/ MiiMi:!, + Mi2Mii \ 
M21M23 +MiiM2i 
M31M33 + M32M34 
\ MiiMi3 + Mi2Mii ) 

Given there exists for each N a set of solutions of ( 5.40| ) such that x' = X- This Floquet set of 
parameters is determined by 

AT 

(5.49) 



where 'im — Z {-M (E) A4) mod(4) is an integer vector. For this set of parameters the propagator Um is such 
that 



N 



-0, 



(5.50) 



where 1^ is the unit operator in [H^] defined in (A. 7). Evidently we cannot change the Floquet parameters 
at each iteration of the map, so (5.4£) defines the quantizability set of the map. For this set of parameters, 
we can then restrict the unitary operators Um to the Hilbert space of the 2L-torus [H^] and thus define 
the corresponding torus propagator as 



N 



(5.51) 



Operators defined through (5.51) inherit the unitary from their plane counterpart and we will call them 
quantum cat maps. 



The fact that the quantum symplectic map commutes with projections on the torus (A. 17) leads to; 

I 



(5.52) 



The unitarity of and ( p).52| ) implies that 



(5.53) 



Having two cat maps M.a and Mb that are quantizable on the same Floquet parameter implies that 

(5.54) 



^ Ma A4f 



^ MaMb' 



that is, the quantization of the composition of different cat maps is equivalent to the composition of the 
quantization of each map. This shows that quantum cat maps form a group, the quantum propagator 
preserving the classical composition laws. This feline quantum group is at the same time a projection onto 
the torus of the metaplectic group and a subgroup of the classical feline group. It then follows that, for any 

power I, the map M'' is also quantizable. 

The quantizability set of the map can be made independent of N if the inhomogeneous term in (5.4C) 
vanishes. This will be the case for matrices M, such that 



m [M^ZM ± a I m =0 mod(2) 



(5.55) 



for any integer vector m. This condition is equivalent to 

L 



'^MijMij+L = even, 



(5.56) 



so that the allowed Floquet parameters are 



X= (M-l)-V. 



(5.57) 
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In conclusion, the Floquet parameters are then denoted by the fixed points of the map (3.61). Tlic QPS can 
thus have any of the fixed points of the map as its origin; without loss of generality, one can take % = 0. For 
the one degree of freedom case ( ^.56 ) implies that the M matrix must be restricted to the family of maps 



even odd 
odd even 



or 



odd even 
even odd 



selected as quantizable by Hannay and Berry Q. 

Let us now provide an explicit construction of quantum cat maps based on the center and chord repre- 
sentations of operators. We start with the quantization of the linear automorphism A4 on the Hilbert space 
[T^e]^ associated to the Euclidean phase space. The linearity of the A4 map implies in the exactness of the 
Van Vleck construction of the propagator IQ . This propagator in the chord representation has the form B 



det 



(5.58) 



where >S'(^) = £"(^,0) is the chord generating function of the automorphism in the plane phase space. In the 
center representation, the Van Vleck propagator is 



Um (x) = 



det 



d^Sjx) 



,iSix) 



(5.59) 



where now S{x) = S{x,0) is the center generating function of the transformation in the plane. Inserting 
( 2.37 ) and ( 2.36 ) for the respective generating functions, we have 



Um{x) = |det(l±aB)r^e*"^" 



and 



Um{0 = |det(l± 5/3)1 



2 p-iiiPi 



(5.60) 



(5.61) 



We may now project these representations on the 2L-torus. For this purpose we must restrict our 
construction to the quantizability set of Floquet parameter x defined through (5.49) and take the chord and 
center symbols of the projected operator. This is performed substituting (5.61) and ( 5.60| ) respectively in 
(A.41) and (A.42). Then we obtain the chord representation of the quantum cat map as an average over 
winding vectors: 



UX,(0-|det(l±3/3)r^(e 



-i2;rAr[i4/3C+i5(/3+0)m+im(/3-3)m+iAm] 



(5.62) 



There are various classical objects present in this formula. Recalling section 2, S{£,, m), the chord generating 
function of the cat map defined in (2.37) appears in the exponent of (5.62). In the amplitude we recognize 
T^, the number of fixed points of the map defined in (2.52), so that. 



(^) = 2^— /e^'27rAr[s(?,m) + .^Am] 



Similarly, for the center representation we obtain 



\Jl^{x) = |dct(l±3B)r^ /t2^N[xBx+xiB-:s)m+im{B+:s)m-^Am] 



(5.63) 



(5.64) 



where we recognize S{x,m.), the center generating function of the cat map (2.36), and r^., the number of 
orbits centered on x defined by ( p. 53 ), so that 



^i2TrN\S{x.m)-^Am\ 



(5.65) 



Both ( |3.63 ) and (5.65) are representations of quantum cat maps of general dimension, showing that 
the quantum propagator is entirely defined in terms of classical objects, except for the term A m in the 
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exponent that denotes the quantum features of the boundary conditions. However, we have to perform 
the average on m to make the representation expHcit. Given that the quantizabihty set (^.4£ ) only admits 
rational values of x, (5.63) and (5.65) are Gaussian sums. The quantizabihty condition (5.55) for the map 
M is equivalent to the condition that the Gaussian sums do not vanish for L = 1. We do not know of a 
similar verification for L > 1. 



In the following we will restrict our attention to maps that fulfill the condition (5.55), such that the 
inhomogeneous term in (5.40) vanishes. As we have already discussed, we can then choose x = without 
loss of generality. The quantum cat maps associated with these Floquet parameters will be denoted by \Jm ■ 
All the examples in section 4 are of this type. 

We now use the periodicity properties ( p. 44 ) of 5(^,111) to show that the exponential function in (5.62) 
is periodic. Indeed, for m' defined in (2.42) and the generating function 5(^,m') defined in (2.44) we have 
that 

^-t2Tr N[S a, m')] ^ g-i2^Ar ,m)] g-i27rAr[C Ak- imFi k- ikAik] ^ {5M) 

where Fi and Ai were respectively defined in ( ^.46 ) and ( 2.47 ). We can now see that g-^^^^li^'^] z 
any integer vector k and for maps that fulfill the condition ( ^.55 ) ^-^'^^^[-^•^^i^-i^'^M — 1^ hence 



1 for 



exp [-i2TrNS (^, m+{M- l)k)] = exp [-i2TrNS{^, m)] , 



(5.67) 



The average in ( ^.63 ) is then periodic, so that, we sum over one period and divide by the number of points 
in the period, to perform such an average. Hence, we must restrict m to the fundamental parallelogram 
defined in (2.50), where there are exactly points m with integer coordinates: 



= 2^{T^y^ E g-»27rAr[iC/3C+iC(/3+0)m+im(/3-3)m] 



(5.' 



meOf 



Since m belongs to O^, the sum in ( ^!68| ) is the sum over the different classical orbits whose chord is ^. 

In a similar way, for the center representation the periodici ty re gion is the parallelogram O^, , defined in 
( 2.51 ), that has exactly points m with integer coordinates ( 2.53|) , so that 



2^T, 



= 2^r, 



E 



J2TTNlS(x,m)] 



^i27rN[xBx+x{B-d)m+^m(B+d)m] 



(5.69) 



i.e. we are taking the sum on the different classical orbits centered in x. The expressions (l^esj) and ds^egj ) 
have exactly the form expected for the semiclassical approximation of the propagator in center and chord 
representations respective ly but in this case they are exact instead of being a mere approximation. 

For generic cat maps, (5.68|) and (5.69) generate Gaussian sums, but in some cases important simplifica- 



tions are possible. Let us start, once more, with the chord representation. We have already seen that the 
matrix f3 has the form f3 = where the barred matrix has integer elements. The simplest case is when 

the cat map is such that the associated /3 matrix itself has integer entries. Recalling that — we will 
transform to an equivalent set of chords, for which the values for ^ are even multiples of r^. So, for any 
chord ^ there is an equivalent chord S, such that 



2t' 



(5.70) 



where the components of S are integer numbers up to N and the components of n are integer up to 2r^ — 1. 
Equation ( 5.70[ ) has solutions with the specified features for any ^, only if N an d 2t^ are coprime numbers. 
We will then restrict N to be an odd integer. In this case, the chords S in ( 5.70|) form a lattice with spacing 

2t' j t 

-j^. A hypercube of side 2t^ has then N x N successive chords S that constitute a basis for translation 
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operators. For the shiiplest case where the matrix f3 has integer entries, the chords S form a lattice of length 
2 with spacing jr. Performing the transformation to chords S, we see that the term g-*2irA'[2=(/3+o)m] _ 



so there is then no S-dependence in the propagator sum (5.6S), leading to(4.6). 
In the same way, the matrix B in the center representation has the form B 

to center points X that are integer multiples of 



Then we transform 



(5.71) 



where the compone nts o f X are integer numbers up to N and the components of j are integer up to 2t^ — 1. 
Again, solutions of ( 5.71 ) with the specified features will exist only if TV and 2t^ are coprime numbers. Thus, 
reflection operators on points X, that form a lattice with separation ^ on a hypercube of side r^, form a 



basis for the Hilbert space of the torus. For center points X in (5.71) , the term g«27rAf[jf (B-3)m] _ -j-j^g 
propagator (5.69) has the form (4.11). 

The above are then special cases where the propagator on the torus has the same form as its equivalent 
on the plane. These cases are then ideal to study quantization, as we have seen in section 4. 

We will now discuss the quantum effect of a similarity transformation of the form 



(5.72) 



As we have already discussed the matrices A4' and A4 represent the same map, but seen on a different 
frame of canonical coordinates. On the torus, we have to restrict both A4 and M to have integer elements 
so that the torus is mapped on itself. This restricts the symplectic similarity transformations to the feline 
transformations. 

The main advantage of the chord and center representation in plane phase space is their symplectic 
invariance[Q. It is well known that linear classical canonical transformations x' — Afx correspond to unitary 
transformations in [TYr] , 



A 



A' = U^rAUJ^\ 



(5.73) 



The effect of such a unitary transformation on the chord and center representation is merely 



A{x)^A{J\fx) and A{£_) ^ A{U£_). 



(5.74) 



Because of the commutation of operator products with projection from the plane to the torus, the effect 

1 ^ 
performed by a quantized cat map on any operator A 



of a similarity transformation A 



TT^ 



that commutes with ijy will be purely classical in the center or the chord representations: 



A{x)^A{JVx) and A{^) ^ A{J\fO- 



(5.75) 



Thus the similarity transformation among quantum cat maps, reduces to the classical similarity transforma- 
tion: 

(5.76) 



so that. 



Vm' {x) = V{J\fx) and Vm' {£.) = U(AAe). 



(5.77) 



However, the quantum cat maps A4' and A4 must be quantized on the same Floquet parameters x- This 
imposes another restriction on the matrix J\f used to perform the similarity transformation, its quantizability 
set must include the Floquet parameter x- That is A4 an d J\f must belong the same quantum feline group. 
For two quantum cat maps U^, and related by ( 5.76 ) and for all power I of the map we have 



Tr 



Tr 



Tr 



(5.78) 
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We have seen in (4.28) that the density of states and hence the spectrum of the system is uniquely determined 
through the traces of the different powers of the map. Then, ( |5.78D shows that U^, and U;J^ related by 



(5.76) have the same quasi-energy spectrum. 

We now show that the quantum period function is also invariant with respect to a feline similarity 
transformation, for thejperiodic case where we can chose the Floquet parameter x = 0. Suppose that the 
QPF corresponding to U^/ and are respectively k^'iN) and kM{N) , then 

(^_M'fMN) ^ ^j^-iMAff-^^^^ =Af-^ {Mf^^^'^J^. (5.79) 

For all points x belonging to the QPS, i.e., a lattice of spacing on the 2L-torus □, we have 

hence 

(M')'=^W = 1 mod(iV), (5.81) 

so that 

kM'{N)<kM{N). (5.82) 
In the same way, for all points x belonging to the QPS 

{Mf^' x=N {M'f^' N~^x = X, (5.83) 

hence 

kM{N)<kM'{N). (5.84) 

Therefore, 

kM'{N)^kM{N), (5.85) 
that is, the QPF of both quantum maps U^/ and coincide. 

6 Conclusions 

In this work we have studied classical and quantum properties of multidimensional Cat maps. In a first step 
the classical study was performed using the symplectically invariant center and chord generating functions. 
They allow us to represent the symplectic matrix by a symmetric one. This is the basis for a complete 
classification of generic four dimensional cat maps. Clearly, the advantage of working with the appropriate 
generating functions will be even more pronounced for cat maps of higher dimension. Two degrees of freedom 
are sufficient to obtain all distinct types of dynamics. Loxodromic behavior appears as a new alternative 
with respect to usual cat maps with one degree of freedom. 

The quantization of cat maps was performed using the recently developed Weyl representation and its 
conjugate chord representation. The semiclassical approximation is exact whatever the number of degrees 
of freedom or the characteristics of the cat map. The spectral properties show the same kind of "pathologies 
" observed for systems with one degree of freedom. Through the quantum periodicity function, we have 
indication of quantal ergodicity and mixing in the semiclassical limit for systems that present this classical 
property. We must note that this is one of the first times that loxodromic behavior is quantized [l5| . 

According to Anosov's theorem, all cases of fully ergodic classical maps are structurally stable, that is, 
a weak nonlinear perturbation leads to a map whose orbits are topologically equivalent to the original cat 
map. The possibility of quantizing such an Anosov map is in no way restricted to one degree of freedom. In 
this way, one can obtain continuous families of quantum torus maps, corresponding to fully chaotic classical 
maps for each type of map ( doubly hyperbolic, loxodromic, etc.). These nonlinear maps will probably avoid 
the spectral anomalies due to quantum periodicity, as was verified for the case of a single degree of freedom 

Acknowledgments: We thanks helpful discussions with J. P. Keating, who provide us Greenman preprint, 
E. Pujals and G. Contreras who provide us Mane's work. We acknowledge financial support from Pronex- 
MCT and A.M.F.R. also thanks support from CLAF-CNPq. This work was partially supported by contracts 
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A The Weyl quantization on the Torus 

In this appendix we develop the mathematical tools needed for the quantization of cat maps, based on our 
previous work on Weyl quantization on the torus H. 



A.l The Hilbert space of the torus 

In a first stage it is important to treat the specification of the Hilbert space of quantum states, or pre- 
quantization, independently from the dynamics of the system. That is, we treat the quantum kinematics, 
corresponding to the geometrical description of phase space at the classical level. Just to simplify the notation 
we limit the presentation for one degree of freedom systems, since the extension for L 7^ 1 is trivial. 
A complete description for prequantization must include Bloch boundary conditions : 

+ = e^'''^''^{q), (A.l) 

*(p+l) e"^"'^''*(p) {A.2) 

where 

= {2TTh)-^/^ J e-''P'''^^{q)dq, (A.3) 

and 2'KiXp and 'i-nixq are fixed arbitrary Floquet angles; that is, the prequantizati on d epends on the vector 
X = (Xpi Xa ) wh ose coordinates are in the range < XaiXa < 1- Solutions to (AJ) and ( A.2 ) with the 
connection (A^) only exist if there is an integer N ^ so that [|l3| 

Then the space of solutions spans a Hilbert space TC^ having the finite dimension N . Two bases of Ti/^ 
appear in analogy with the position and momentum eigenvectors, here denoted by, |q„ > and |pm > with 
n,m = 0, 1, . . . , — 1. According to the Bloch type boundary conditions, we will then define the position 
state on the torus as an average over equivalent positions in the plane phase space: 

Xq 

with the Hermitian structure 



qn>=(l^^ + A:>e^-'=>^0 , (A.5) 



< q™|q„ >N,x=< qm|q« <5We'^('"-")'^^ (A.6) 



The identity operator in Ti/^ is obtained through 

Af-l 



1^ = 1^™ >< (^■'') 

Tn=0 

In analogy, the momentum eigenvectors 

Ip™ >= (1^^^ + k > e-^-^^^^'j^ , (A.8) 

normalized such that 

< P™|p„ >= (5We-^(™-")x,. (A.9) 
The bases are exchanged with the transformation kernel, 

forming a unitary matrix (finite Fourier transformation). According to ( [A.4| ) the classical limit corresponds 
to A^ ^ 00 . 
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We then see that positions and momenta form a discrete web on the torus, that we wih call from now 
on the quantum phase space QPS in accordance with For the case of L degrees of freedom, the Hilbert 
space of the 2L-torus is [H%] ='H%xn^ X ... X H]^ [L times). For each degree of freedom there is a grid 



structure, as above, so the dimension of the Hilbert space 

[T-l%] is N^- In this case the Floquet parameter 



is a 2_L-dimensional vector. 



A. 2 Projector of the plane on the Torus: 

We envisage the existence of two Hilbert spaces, one associated to the torus, 7i^, and another, TYr, associated 
to the Euclidean phase space, formed by the square-integrable distributions that we extend to quasiperiodic 
distributions. It is easy to see that Ti^ is the subspace of TYr that obeys the boundary conditions (A.l).We 
then define a projection of TYr on through the operator 1^, 



where 



|*>=i^|V>, 

N-l 

1^ = X! i*^" X 



(A.ll) 
(A.12) 



n=0 



We see that 1^ is Hermitian 1^ 



1^ and, inserting ( A. 6 ), we obtain that it is a projector 



J- AT- 



(A.13) 



For all operators A acting in TYr there is a projected torus operator which acts on TC^ 



AX 



1^ 



(A.14) 



Let us now suppose that the operator A transforms TL^ into itself, i.e. 



then 



(A.15) 
(A.16) 



and A is said to be a torus invariant operator. 
For two operators of this kind: 



J- AT ^ J- AT -O J- AT 



(A.17) 



So, the product of the projected operators is the same as the projection of the product of operators. For many 



degrees of freedom, we will simplify the notation by using 1^ to denote both the projector on [7i^] 



and 



Tiff. No confusion is made since we will always specify the number L of degrees of freedom we are working 



A. 3 Restriction of the translations and reflections to the torus 

Projecting the translation and reflection operators on the torus we obtain M 



rpX 





if there are r and s integers so that f = iw^w) 
otherwise 



(A.18) 



where the torus translation operators T|- = T^^ are defined through 



(A.19) 
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Their interpretation as Translation operators in QPS is clear with this last expression. The x dependence is 
from now on implicit. 

Performing m = ( ''1 integer loops around the irreducible circuits of the torus we get 

\ 'm-q J 

T„ = e*2^^[(*)^'»+i™5m] ^ (A.20) 
so that the symmetries of the translation operator are 

T^+rn = e"''^[('"^)^'"+^"^'"]Tj. (A.21) 

To obtain a basis of operators, r and s must then run in [0,N — 1], that is, we only need translations that 
perform less than one loop around the torus. 

Similarly, projecting reflection operators, we have 



1^ 7? 1^ 

-■-JV a:-'-JV 



Hxa.b if there are a and b semi-integers so that x = ^^^p-, ^7^) 22) 



otherwise 



where 

2N-12N-1 

^-.^ = ^ E E e-^'-^^-^% (A.23) 
arc the torus reflection operators on the center point Xa,b- The unitarity of Ra, is ensured by the action of 



the operators on the Hilbert space 



Rx|q„ >= e^^^^^-^^^^+^'-^lqsb-n > and R^|p„ >= e'^2(a-m)(6+x<,)|p2^_^ > . (A.24) 

Then R^, reflects the QPS web about the point x = ( , ) = {xp,Xq). We then need to include 
half-intcgcr values of a and b to perform these reflections. Symmetry properties for reflection operators on 
points differing in half loops are, 

R^_^» = (■_l)b'»p+am,+mpm,JV^_^ ^ gi27riv[^ Am+im5m] ^^ 25) 

The N"^ independent operators needed for a basis of are obtained with (a, b) belonging to 0, 1/2, ... , ^^^^,i.e., 
only one quarter of the torus is needed. This is named the Weyl Phase Space WPS. The traces of the trans- 
lation and reflection operators are: 

Tr(Tj) = A/-e'T^('r+'-Xp-«Xp)^Wj(^) = n e*T^(¥+'^x,-%) (a.26) 

and 

TV(R,) = /^(x) = i(l + (-l)2- + (-l)2& + (-l)2-+26+A') (A.27) 



if N is even and a or b half-integers 
2 if N is even and a and b integers 

1 if N is odd and a or 6 integers 
—1 if N is odd and a and b half- integers 



(A.28) 



A. 4 Operators and their Symbols 

With the set of torus translation and reflection operators, we can represent any operator in . The chord 
representation of an operator A is defined as; 

A(0 = Tr (at . (A.29) 



31 



From the symbol, we recover the operator 

A-]^E Am^-^A(OT,. (A.30) 

r,s=0 { 

Performing m foops around the torus, the symbol becomes 

A(^ + m) = e'2^^[(i-*)^"^+3™5m] A(0. (A.31) 
The center representation on the torus of an operator A is defined as 

A(x) = Tr (^AR^rj . (A.32) 
From the symbol we recover the operator through: 

N-l 

^ = ^ E ^-A(^) ^ ^ E ^^M^)- (A.33) 

a,fc— X 



The symmetry properties of R^, (A. 25) implies in , 

From the chord or center representations we obtain the more familiar position representation of the 
propagator as 



A(q,„,q„) = - E A(e„6n-n)e^'^^«-('^\ (A.35) 



A(q„,q„) = ^ E A(xp,x^)e^2-^-''(«'"-«"). (A.36) 



N 

x„=0 



The trace of operators in TC^ are obtained as 



N 

N- 



Tr 



(a) = A(e = 0) = 5] A{x)Mx) = IJ2 



a,fc=0 



The representation of the identity operator on the torus Hilbert space Ti.^ has the form 

l{x) = fN{x) and li^) = NSf\ (A.38) 
For the chord representation to denote a unitary operator U, it has to fulfill the condition 



while the center representation of unitary operators is restricted to symbols such that, 

2 



(A.39) 



(^) E U(xi)U*(x2)e'4"^("""i)^("-"i) = l{x) = fNix). 



(A.40) 
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A. 5 Relation between symbols 

We shall now describe how the symbols on the torus can be obtained from their counterparts on the plane. 
The result is directly given for L degrees of freedom. Starting with the chord representation, the torus 
symbol at points = where ^ is a 2L integer components vector, is calculated in so that 

\ / m 

Note that we have to perform a phase weighted average on equivalent points to obtain the symbol on the 
torus. In a similar way the symbols in the center representation on points x — ix + x) E^re also calculated 
in j|] resulting in 
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